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For  the  past  decade,  extensive  mathematical  modelling  has  been  conducted  on  the  design  and  optimiza¬ 
tion  of  liquid-feed  direct  methanol  fuel  cells  (DMFCs).  Detailed  modelling  of  DMFC  operations  reveals 
that  a  two-phase  flow  phenomenon  at  the  anode  and  under-rib  convection  due  to  the  pressure  dif¬ 
ference  between  the  adjacent  channels  both  contribute  significantly  to  mass-transfer  in  a  DMFC  and 
its  output  performance.  In  practice,  comprehensive  simulations  based  on  the  finite  volume  technique 
for  two-phase  flow  require  a  high  level  of  numerical  complexity  in  computation.  This  study  presents  a 
complexity-reduced  mathematical  model  that  is  developed  to  cover  both  phenomena  for  a  realistic,  but 
fast,  in  computation  for  the  prediction  and  analysis  of  a  DMFC  prototype  design.  The  simulation  results  are 
validated  against  experimental  data  with  good  agreement.  Analysis  of  the  DMFC  mass-transfer  is  made 
to  investigate  methanol  distribution  at  anode  and  its  crossover  through  the  proton-exchange  membrane. 
From  a  comparison  of  the  influence  of  two-phase  flow  and  under- rib  mass-transfer  on  DMFC  performance, 
the  significance  of  gas-phase  methanol  transport  is  established.  Simulation  results  suggest  that  both 
the  optimization  of  the  flow-field  structure  and  the  fuel  cell  operating  parameters  (flow  rate,  methanol 
concentration  and  operating  temperature)  are  important  factors  for  competitive  DMFC  performance 
output. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  considered  to  be  a  promising-next  generation  power  source,  especially  for  portable  and  micro¬ 
scale  applications,  due  to  its  simplicity,  low  temperature  operation,  and  relatively  high  specific  energy  [1,2].  Nevertheless,  the  present 
performance  of  DMFCs  is  interior  to  that  of  state-of-the-art  lithium-ion  batteries.  The  key  barriers  to  its  rapid  commercialization  are 
poor  methanol  oxidation  kinetics  and  methanol  crossover  through  the  polymer  membrane  [3,4].  Although  extensive  research  is  being 
conducted  world-wide  on  new  materials  for  performance  improvement,  other  work  on  DMFC  design  and  optimization  using  currently 
available  materials  emerge  is  also  an  important  topic  [5].  The  availability  of  a  mathematical  model  of  a  DMFC  helps  to  understand  the 
interactions  inside  the  cell.  Furthermore,  a  model  serves  as  a  tool  to  optimize  DMFC  design  and  allows  more  accurate  prediction  of  the 
performance  before  a  prototype  is  made  [6].  Therefore,  a  good  mathematical  model  is  essential  for  the  design  and  development  of  a 
competitive  DMFC. 

In  the  development  of  DMFC  models,  the  mass  transport  of  reactants  has  been  found  to  be  an  important  and  crucial  issue.  Accordingly, 
their  parameter  has  attracted  more  attention  in  model  development  [7,8].  In  early  work,  only  single-phase  mass  transport  in  one  or  two 
dimensions  was  adopted  in  models  [9-13].  The  typical  assumption,  such  as  specified  in  the  models  by  Meyers  and  Newman  [9-11]  and 
Baxter  et  al.  [12],  was  that  carbon  dioxide  produced  in  the  reaction  remains  fully  dissolved  in  the  liquid  fuel.  This  assumption  was  found 
to  be  inadequate  by  studies  of  the  flow  visualization  of  carbon  dioxide  growth  and  its  ejection  from  the  diffusion  layer  [14-17].  Therefore, 
unlike  other  types  of  fuel  cell,  methanol  transport  in  a  DMFC  occurs  in  a  liquid-gas  two-phase  flow,  which  consists  of  aqueous  methanol  and 
gaseous  carbon  dioxide.  Accordingly,  the  mass  transport  process  is  more  complex  and  it  makes  the  development  of  a  two-phase  model  rather 
challenging.  Wang  and  Wang  [18-20]  successfully  simulated  the  two-phase  flow  and  transport  in  a  DMFC  by  applying  their  multiphase 
mixture  (M1 2)  flow  model.  Capillary  action  was  believed  to  be  the  dominant  mechanism  for  two-phase  flow  within  the  diffusion  layer.  The 
simulation  clearly  showed  the  effect  of  the  gas  phase  transport  of  methanol  on  the  DMFC  output  performance.  Thus,  two-phase  flow  is 
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Nomenclature 

Symbols 

C  molar  concentration  (mol  cm-3 ) 

D  diffusion  coefficient  (cm2  s-1 ) 

dh  hydraulic  diameter  (cm) 

F  Faraday  constant  (96,487  As  mol-1) 

hm  mass-transfer  coefficient  between  channel  and  diffusion  layer  (cm  s-1 ) 
I cen  current  density  (A  cm-2 ) 

j  current  density  per  unit  volume  (A  cm-3 ) 

I<  proportional  coefficient  for  methanol  electro-osmosis  drag  coefficient 

k  permeability  of  porous  materials  (cm2 ) 

kH  Henry’s  law  constant  (Pa) 

l  channel  length  (cm) 

M  molecular  weight  (g  mol-1) 

N  mass  flux  (mol  (cm2  s)-1 ) 

nEOD  water  electro-osmosis  drag  coefficient  (cm3  mol-1 ) 

Piniet  inlet  flow  pressure  (Pa) 

P  pressure  (Pa) 

A P  pressure  drop  (Pa) 

R  universal  gas  constant  (8.314J  mol-1  K-1 ) 

Rcontact  ohmic  contact  resistance  ( £2  cm2 ) 

S  liquid  saturation 

Sh  Sherwood  number 

T  temperature  (I<) 

Vceu  cell  voltage  (V) 

v  flow  velocity  (cm  s-1) 

x  coordinate  along  channel  (cm) 

y  coordinate  across  channel  to  membrane  (cm) 

Superscripts 
MeOH  methanol 

C02  carbon  dioxide 

H20  water 

02  oxygen 

Subscripts 
A  anode 

C  cathode 

CH  channel 

CH-DL  channel  to  diffusion  layer 

Cov  convection 

DL  diffusion  layer 

Diff  diffusion 

EOD  electro-osmosis  drag 

eff  effective  value 

g  gas  phase 

in  inlet 

l  liquid  phase 

m  membrane 

ref  reference  value 

sat  saturated 

U-rib  under-rib 

V  vapor 

Greek  letters 

aA  anodic  transfer  coefficient  at  anode 

ac  cathodic  transfer  coefficient  at  cathode 

a  gas  void  fraction  in  channel 

1 7  overpotential  (V) 

s  porosity 

r  tortuosity  factor 

Km  proton  conductivity  of  PEM  (S  cm-1 ) 
a  interfacial  tension  (S  cm-1 ) 

li a  viscosity  (kg  (cm  s)-1 ) 
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p  density  (gem-3) 

X^qdH  methanol  electro-osmosis  drag  coefficient  (cm3  mol-1 ) 
^MeOH  liquid  methanol  molar  fraction 


recognized  to  be  a  more  realistic  model  for  mass  transport  in  DMFCs  [21  -25].  The  effect  on  mass-transfer  by  under-rib  convection  within  the 
diffusion  layer  was  also  investigated  extensively,  although  mostly  in  hydrogen  proton-exchange  membrane  fuel  cells  (H2  PEMFCs)  [26-31  ]. 
Under-rib  convection  is  induced  by  the  pressure  drop  between  two  adjacent  channels.  It  is  understood  that  such  convection  dominates 
lateral  species  transport  in  the  diffusion  layer.  Recently,  under-rib  mass  transport  was  found  to  have  a  substantial  influence  on  DMFC 
performance  [7,32]. 

Most  of  the  available  mathematical  models  are  applied  numerically  using  the  finite  volume  technique,  together  with  the  computable 
fluidic  dynamics  (CFD)  method.  Detailed  analysis  on  the  species  transport  is  made  by  the  complete  Navier-Stokes  equation  set,  but  involves 
vast  numerical  complexity. 

In  this  study,  a  concise  mathematical  model  is  developed  for  quick  numerical  simulation  to  allow  performance  prediction  and  analysis 
of  a  DMFC.  Two-phase  flow  at  the  anode  and  pressure-driven  convection  are  both  included  in  the  proposed  model  and  thereby  enable 
accurate  and  realistic  performance  prediction  and  analysis.  The  results  are  validated  by  laboratory  measurement  data. 


2.  Mathematical  model 

A  schematic  of  a  DMFC  with  a  single  serpentine  channel  (SSC)  is  shown  in  Fig.  1.  The  anode  and  cathode  sides  have  a  similar  structure 
that  consists  of  a  fuel  channel  and  diffusion  layer.  A  Nation®  type  proton-exchange  membrane  (PEM)  is  sandwiched  between  the  anode 
and  the  cathode,  with  catalysts  supported  on  carbon  particles  on  each  side  for  the  fuel  cell  reactions,  namely  methanol  oxidation  on  a 
platinum-ruthenium  (Pt-Ru)  alloy  at  the  anode  and  oxygen  reduction  on  a  platinum  (Pt)  at  the  cathode.  The  electrocatalytic  reaction  at 
the  anode  can  be  expressed  as 

CH3OH  +  H20  -*  C02  +  6H+  +  6e-  (1) 


Fig.  1.  Schematic  structure  of  DMFC  with  SSC  pattern  for  both  sides.  Face  A  represents  a  cross-sectional  view  of  region  with  under-rib  mass  transport.  Face  B  represents  the 
region  without  under-rib  mass  transport. 
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and  the  reaction  at  cathode  is 

6H+  +  6e“  +  1 .502  -*  3H20  (2) 

For  simplicity,  both  catalyst  layers  are  assumed  to  be  infinitely  thin  at  the  interface  between  the  diffusion  layer  and  the  PEM  in  the 
proposed  model.  Accordingly,  the  methanol/oxygen  diffusion  with  electrochemical  reaction  is  assumed  ideal  in  these  layers.  The  modelling 
work  on  cell  performance  and  mass  transport  therefore  addresses  the  flow  channel,  diffusion  layer  and  PEM  regions. 

The  SSC  is  treated  in  two  regions  for  mass  transport  analysis,  as  shown  in  the  top-view  of  the  entire  flow-field  in  Fig.  1.  The  under¬ 
rib  convection  mass  flux  takes  place  between  neighboring  channels.  In  this  region,  species  flux  is  not  only  driven  by  the  concentration 
gradient,  but  also  by  the  pressure  difference  between  channels  according  to  Darcy’s  convection  flow.  The  liquid  pressure  in  the  cross- 
section  is  assumed  to  be  homogeneous  in  the  channel  and  beneath  diffusion  layer.  For  the  region  with  a  horizontal  short  turning  of  flow 
channel,  the  mass  flux  is  assumed  to  be  driven  mainly  by  the  concentration  gradient  and  the  convection  flux  can  be  neglected  due  to  the 
relatively  small  pressure  drop  along  the  channel.  The  phenomenon  is  indicated  in  the  cross-sectional  view  of  face  B  in  Fig.  1.  The  two-phase 
mixture  flow  is  the  main  transport  mode  for  methanol  and  carbon  dioxide  at  anode.  In  the  anode  channel  region,  the  laminar  two-phase 
flow  is  treated  homogeneously  and  similar  gas  and  liquid  phase  velocities  are  maintained.  At  the  cathode,  the  two-phase  effect  is  not 
considered  due  to  the  very  small  amount  of  liquid  water  [18].  The  crossover  of  methanol  through  the  PEM  from  the  anode  to  the  cathode 
is  assumed  to  be  oxidized  instantaneously  at  the  cathode  catalyst  layer  [33,34]. 

Other  assumptions  for  the  mass  transport  in  the  model  development  are:  (1 )  all  processes  in  the  DMFC  system  are  isothermal  and  under 
a  steady-state  condition;  (2)  the  liquid  and  gas  considered  in  the  model  are  incompressible  and  their  flows  are  laminar  due  to  the  low 
Reynolds  number;  (3)  porous  diffusion  layers  are  isotropic  and  homogeneous  and  characterized  by  effective  porosity  and  permeability. 

2.1.  Methanol  crossover  in  PEM 


In  a  PEM,  the  crossover  methanol  flux  arises  in  the  presence  of  excessive  un-reacted  methanol  at  the  diffusion  layer]  PEM  interface  due  to 
slow  methanol  oxidation  kinetics.  The  crossover  is  driven  by  diffusion  under  a  concentration  gradient,  electro-osmosis  drag  proportional  to 
cell  current  density  and  local  methanol  concentration,  as  well  as  a  pressure  difference  and  local  methanol  concentration  induced  convection 
across  the  PEM  [35],  i.e., 


NMeOH  =  _Dl 


rdC^°» 


™’eff  dyn 


,  vMeOHrA/ 
A  EOD  ^  rri 


i  1ceU  _  i  p  HV 

p 


r  APrr 


(3) 


where  the  methanol  electro-osmosis  drag  coefficient  is  x^odH  =  ^o^n£OD  anc*  effective  diffusion  coefficient  is  =  (sm)Tm D^e0H t 

respectively. 

The  flux  is  constant  in  the  PEM  at  steady  state  and  therefore  the  methanol  concentration  distribution  can  be  solved  with  the  boundary 
conditions  of  Cj%e0H  =  \y_t  at  the  interface  between  anode  catalyst  layer  and  PEM,  and  Cj jfe0H  =  0  at  the  cathode  catalyst  layer 

assuming  that  all  crossover  methanol  molecules  are  consumed  by  combination  with  oxygen  at  cathode  side;  i.e., 

C„»H  c„„„,  “P«r»'/i>S3>)y»)-«xp((r;/DS»»)t„)  APm 

"  " U-m  i -exp((r;/D“<“»„)  •'«-*«»  f  tm  m 

Therefore,  the  crossover  methanol  flux  in  PEM  is  given  as 

exp((r;/D^H)tm) 


N! 


r  MeOH 


=  c; 


■MeOH  I 


A—DL  |  y=tfy_ 


-dl  Aexp ((r;/D^H)tm)_1 
Accordingly,  the  diffusion,  electro-osmosis  drag,  as  well  as  convection  partial  flux  are 


(5) 


i vfMeOH  _  _  rMeOH 
iym,Diff  ~  nr 


exp((rAVP^)tm) 

A  i  d  /  nMeOH  \ 


A~DL  A  1  -  exp((r^/D^)t, 


NMeOH  __  vMeOH  Leg  r MeOH 
J  mn  —  A  EOD  ^ A  nl 


PjMeOH  =  _pP  ^rn^MeOH 


A  /  m,eff 
lA* 


exp  ((r^/D^OH 


)tm )  -  exp((  )tm ) 


y=tA_DL  1-exp  ((rJ/D^ptn 


(6) 


m,Cov 


^ A-DL 


exp  ((r^/Df  o«)ym)  -  exp((  rj  /  )tm ) 


y=tA-DL 

2.2.  Methanol  mass  transport  at  anode 


m,eff _ 

l-exp  ((rJ/D%™)t„ 


At  the  anode,  methanol  solution  is  supplied  from  the  anode  flow  channel.  It  is  treated  differently  with  and  without  under-rib  convection 
in  the  corresponding  region,  as  follows. 


2.2.1  Under-rib  mass  transport  region 

In  the  region  with  under-rib  mass  transport,  the  Darcy-type  methanol  flow  consists  of  diffusion  flux  and  a  convection  under-rib  flux 
governed  by  the  Darcy  model,  i.e., 

rfpMeOH  \,P  A  D 

\jMeOH  _  nMeOH  A-DL  KA-DL  rMeOH  1X1  U-rib 

NA-DL  -  UA-DL,eff  dyA  DL  ^  ^A-DL  ^ 
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The  effective  diffusion  coefficient  of  methanol  D^eOH^  in  the  layer  also  relates  to  the  porosity  and  tortuosity  of  the  porous  region  similar 
to  that  in  the  PEM.  Mixed  parameters,  which  will  be  elaborated  in  the  two-phase  section  below,  will  be  applied  assuming  that  methanol  is 
transported  in  both  the  gas  and  the  liquid  phases  in  the  diffusion  layer. 

The  pressure  drop  APu_rib  in  the  above  flux  expression  is  the  difference  in  the  local  pressure  at  corresponding  points  between  adjacent 
channels.  Assuming  laminar  flow  in  the  fuel  channel,  this  pressure  drop  of  methanol  flow  along  channel  can  be  derived  from  laminar  flow 
theory  [36],  i.e., 

A  r.-c„(x)=32^-f^x  (8) 

dh 

Therefore,  the  pressure  difference  between  neighboring  channels  is: 


A Pu-ribW  =  [(Pfntet  -  PM)  -  (. Pinlet  ~  PM))]  =  ^zCMA[2 n/„  +  (2n  -  1)1*  -  2x] 


dl 


where  n  is  the  parallel  channel  number. 


Accordingly  : 


ArMeOH 

\rMeOH  _  nMeOH  A-DL  ,  r2  rMeOH 

nr  -  uA  DLeff  dyA  DL  ^'a^a-dl 


M-DL 


=  32uA  CHkpA_DL[2n^  +  (2n  _  1  )h  _  2x] 
lhdi 


(9) 


(10) 


Following  mass  conservation  /dyA-DL  =  0,  the  distribution  of  methanol  concentration  across  diffusion  layer  is  given  by 


rMeOH 
^ A-DL 


_  rMeOH 
y=0  A-DL 


exp  )yA_DL)+  CMeOH 

y=tA-DL  J 


y=tA-DL 


-  CA-DH1  L0exp((/7/Df ; ™e// W) 


1-exp  ((^/D^OHeff)tADL) 

- A-DL  ~  \ A-DL  \y=tA_DL  at  the  an0de  Cata'VSt  laVer  and  CA-DL  =  CA-DL  |y=0 

anode  channel  and  diffusion  layer.  The  methanol  mass  flux  at  anode  diffusion  layer  is  then  described  as 


The  boundary  conditions  are 
and  diffusii 

p2  rMeOH  I 

jala-dl  |y=tA DL 


(11) 


at  the  interface  between  the 


AjMeOH 

nA-DL 


~  n  W  |y=0exp ((r,2/Df  °«efft,Di)) 


1-exp  ((r2/DM°»efftA_DL)) 


=  ^+<e0H 


(12) 


which  corresponds  to  the  electricity  generation  as  well  as  the  methanol  crossover  through  the  PEM. 

Hereafter,  a  first-order  electrochemical  reaction  is  adopted  in  the  proposed  model  for  both  methanol  oxidation  and  oxygen  reduction 
hereafter.  Consequently,  Tafel  kinetics  of  the  first  order  are  employed  to  describe  the  reaction  current  of  methanol  oxidation  at  the  anode 
catalyst  layer,  namely: 


,•  _  r3  rMeOH 
Ja-  *  A  s 


a  A-DL  \y=tA_DL 


FI  =  U^0H/C^0H)exp((aAF/RT)r]A) 


(13) 


Combining  Eqs.  (5)  and  (11),  the  equilibrium  relation  between  methanol  concentration  at  catalyst  layer  and  channel|diffusion  layer 
interface  can  be  found,  i.e.: 


rMeOH  I  _  a/W,  U-rib  rMeOH  I 

A-DL  I  y=tA-DL  A-DL  A-DL  |y=Q 


(14) 


where  the  equilibrium  constant  's  defined  as 


ajW,  U-rib 
* A-DL 


n  eXP( ( Ct  /^A-DL^ff  )ta-DL )/( 1  -  exp ((^/D^OHeff)tA DL)) 


MU  1  -  exp((r,2/Df  ™eff)tA-DL))  +  M}  exp((r;/D^)tm)/l  -  exp((r;/D^)tm))  -  (rA3/6F) 


(15) 


Due  to  the  continuity  of  methanol  flux  at  the  anode,  the  flux  in  the  diffusion  layer  can  be  described  by  effective  mass-transfer  at  its  interface 
toward  the  anode  flow  channel: 


AjMeOH  _  juMeOH 
nA-DL  ~  nA-CH-DL 


_  uMeOHf  rMeOH 
~  m,eff  ^ A-CH 


rMeOH 
L A-DL 


ly=0 


) 


(16) 


where  is  the  effective  mass-transfer  coefficient.  As  an  approximation,  the  effective  diffusion  coefficient  relating  to  a  two-phase  effect 
of  fully  developed  laminar  flow  can  be  expressed  as  [18] 


r\MeOH 

hMeOH=Sh  A~CH’eff 

m,6ff  tA-CH 


(17) 


As  a  result,  the  equilibrium  relation  between  methanol  concentration  at  the  anode  diffusion  layer  and  its  channel  interface  can  then  be 
derived  by  combining  Eq.  (12),  Eq.  (14),  and  Eq.  (16): 


rMeOH  I  _  AyW, U-rib  rMeOH 
^ A-DL  |y_o  -  ^A-CH-DL^A-CH 


(18) 
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where  'PYch-dl  *s  an  equilibrium  constant  under  the  assumption  of  instantaneous  equilibrium  at  the  interface,  i.e.: 


9) 


w,  U-rib 


uMeOH 

lm,eff 


DL  /?(*£{£*  -  exp^/Df °D»eff)tA-DL)n  -  exp ((^/Df ^e//)tA_DL))  + 

According  to  mass  conversation,  the  methanol  flux  from  the  channel  to  the  diffusion  layer  can  be  described  as 


-tA-CHVA-CH 


dx 


_  i \jMeOH  _  uMeOH  a  ,/vW,  U-rib  \^MeOH 

~  A-CH-DL  ~ 


m,eff 


M  _  \T/W,u-rw  \rivieUH 
^  1  * A-CH-DL^A-CH 


(19) 


(20) 


The  concentration  distribution  of  methanol  along  the  anode  channel  Cjj^f(x)  can  be  obtained  by  solving  the  differential  equation  Eq 

y=tA-DL 


(20).  Subsequently,  its  distribution  at  the  channel] diffusion  layer  interface  | 


(x)  and  the  anode  catalyst  layer  |  _  (x)  can 


be  derived. 


2.2.2.  No  under-rib  mass  transport  region 

Similar  to  that  of  the  under-rib  mass  transport  region,  the  mass  flux  of  methanol  is  mainly  driven  by  diffusion  and  the  under-rib 
convection  is  neglected: 


f]r'MeOH 

AjMeOH  _  nMeOH  uy~A-DL 

A~DL  ~  U^eff  dyA_DL 


(21) 


Based  on  mass  conservation  dN^f™ /dyA_DL  =  0,  the  methanol  concentration  is  found  linearly  distributed  in  the  diffusion  layer: 


_  rMeOH 
" A-DL  I  y=tA_DL  LA-DL 

tA-DL 


y-°7  ..  ,  rMeOH 

yA-DL  +  tA-DL 


y= o 


Similarly,  the  equilibrium  relations  of  methanol  concentration  at  different  locations  are  given  as 


w/o, U-rib 


-A-DL  | y=tA_DL  -  * A-DL 


rMeOH  I  _  tf/W/o, U-rib  rMeOH 
V A-DL  |y=o  ~~  A-CH-DL  ^ A-CH 

where  the  equilibrium  constants  are 


y=0 


(22) 


(23) 


^w/o, U-rib 
^ A-DL 


^w/o, U-rib 
^ A-CH-DL  : 


(  ^A-DL,  eff  /  ^A-DL  ) 


r;exp((r;/D^H)tm 


_( riMeOH  /tii  ~  ^  m,eff 

A-DL,eff/  A-DL)  ^  _  exp((/^l /DMeOH)tj 


r2 

1  A 


,)  6F 


(24) 


m.eff 


(i  -  /«-»> + 

As  a  consequence  of  mass  conservation  of  methanol  in  the  channel: 
dC™°" 


-tA-CHVA-CH  - 


A-CH 

dx 


_  uMeOH _  n/W/o, U-rib, rMe 
~  lm,eff  '  1  ^ A-CH-DL  '^A- 1 


w/o, U-rib,  rMeOH 


A-CH-DL  J^A-CH 


(25) 


As  a-ch-dl  not  related  to  the  methanol  concentration  in  the  channel  Cj^SP,  the  final  solution  for  methanol  distribution  Cj^2?(x)  in 


A-CH-DL 

the  channel  is 


Cfg?(x)  =  Kexp(- 


uMeOH (a  _  y^w/o, U-rib, 
,  m,eff  ^  1  ^ A-CH-DL  ' 


tA-CHVA-CH 


X) 


(26) 


where  I<  is  the  constant  determined  by  the  boundary  conditions  for  the  region  without  under-rib  mass  transport  in  the  entire  flow-filed. 
Subsequently,  the  distribution  of  methanol  in  the  diffusion  layer,  |^_0(x) anc*  |y_t  can  derived. 


The  cell  current  density  is  calculated  by  integrating  the  local  current  density  of  Eq.  (13): 

/t total 

C?«|y=t4DWdx 


(27) 


where  C^e^\  (x)  is  previously  derived  for  regions  with  and  under-rib  mass  transport. 

n  UL  'y=tA-DL 

The  effective  methanol  diffusion  coefficient  DA-DLeff  and  mass_transfer  coefficient  h^eff  remain  undetermined  as  these  two 
parameters  are  affected  by  two-phase  transport  due  to  the  existence  of  carbon  dioxide. 


2.3.  Two-phase  effect  at  anode 

The  diffusion  layer  used  in  a  DMFC  is  normally  made  of  carbon  paper  with  a  microfibrous  structure  and  a  pore  size  that  ranges  from  a  few 
to  tens  of  microns.  The  paper  is  usually  teflonized  to  impart  a  hydrophobic  property.  Together  with  carbon  dioxide  generation,  liquid-gas 
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coexist  in  the  capillaries  of  the  layer.  The  two-phase  transport  in  the  diffusion  layer  is  then  governed  primarily  by  capillary  force  due  to 
the  surface  tension  of  the  liquid-gas  interfacial  meniscus  in  the  porous  structure  [37].  Therefore,  the  contact  angle  and  liquid  retention  are 
important  parameters  in  this  region.  The  capillary  pressure  between  the  gas  and  liquid  regions  can  be  typically  expressed  as 


Pc  =  pg  Pi  =  a  cos  0 


1/2 

F(S) 


(28) 


where  6  is  surface  contact  angle,  and  the  Leverett  function  F(S)  is  the  dimensionless  capillary  pressure  as  a  function  of  the  liquid  saturation, 
S,  which  describes  the  ratio  of  liquid  volume  within  the  porous  structure,  namely: 

F(S)  =  [1 .417(1  -S)- 2.120(1  -  S)2  +  1.263(1  -S)3]  (29) 


As  the  flow  is  assumed  to  be  homogeneous,  the  liquid  saturation  S  takes  the  average  value  of  the  entire  flow  region  for  the  proposed 
model. 

With  back  transfer  of  carbon  dioxide  through  the  diffusion  layer  to  the  anode  channel,  a  localized  thermodynamic  equilibrium  occurs 
in  the  liquid-gas  region.  Therefore,  the  gas  phase  at  the  anode  diffusion  layer  also  consists  of  saturated  water  and  methanol  vapor,  i.e.: 


n  _  pH20  ,  pMeOH  ,  nC02 
1  S  ~  1  V,sat  +  1  V,sat  +  1  g 


(30) 


As  a  result,  water  and  methanol  are  treated  as  a  two-phase  mixture.  Hence,  following  two-phase  mixture  parameters  will  apply. 


Density  :  p  =  piS  +  pg(  1  -  5) 


Concentration  :  C  =  C/S  +  Cg(l  -  S) 


(31) 


Diffusion-coefficient :  pD  =  PiSDt  +  pg{  1  -  S)Dg 

The  effective  diffusion  coefficient  for  each  phase  of  methanol  can  be  expressed  as 

da-dl,i  =  (sAS)TAD^e0H  and  D“fg£g  =  [eA(l  -  S)]TAD^e0H  (32) 

The  water  saturation  vapor  pressure  Pyjg  can  be  readily  found  from  a  steam  table  [38],  whereas  Henry’s  law  is  employed  for  the 
methanol  vapor  pressure  [18]. 

PV,sat  ~  kHXi  (33) 

where  x^e0H  is  the  liquid  phase  methanol  molar  fraction.  For  an  average  diluted  solution,  this  fraction  can  be  determined  by  the  average 
liquid  methanol  concentration  in  the  diffusion  layer,  i.e.: 


„MeOH 

A/ 


rMeOH 
^ A-DL,l 


CH2U  >>cMeOH  rMeOH 
A-DL,l~A~DL,l  LA-DL,l 


rH20  ,  rMeOH 
A-DL,l  ^  ^ A-DLf 


CH  2^ 

^ A-DL,l 


A  aH20  rMeOH 
1V1l  ^ A-DL,l 

P^° 


(34) 


By  applying  the  ideal  gas  law  as  well  as  the  two-phase  mixture  parameter  definition,  the  methanol  vapor  pressure  can  be  related  to  the 
average  methanol  concentration,  as  follows: 

pMeOH  _ /<h  •  M^2° / p^2° _ cMeOH  fon 

v,sat  5  +  [/<HMf20(l  _  S)/(p^2°RT)]  A~DL  1  j 

The  carbon  dioxide  is  assumed  to  be  non-soluble  in  liquid  to  make  the  model  simpler.  Therefore,  its  partial  pressure  can  be  related  to 
the  average  concentration  Cj^L  in  diffusion  layer  by  employing  ideal  gas  law,  P£° 2  =  RTC '^%L. 

In  the  channel  region,  two-phase  transport  is  rather  simple  as  a  homogeneous  bubble  flow  is  assumed.  The  gas  phase  carbon  dioxide 
is  saturated  with  water  and  methanol  vapor.  The  total  pressure  of  the  gas  phase  species  is  equal  to  that  of  liquid  due  to  an  assumption  of 
homogeneous  steady  flow.  Similarly,  the  partial  pressure  of  each  species  is  proportioned  to  its  average  concentration  in  the  anode  channel, 
i.e.,  and  C^2H  respectively.  Assuming  a  is  a  void  fraction  which  represents  the  gas  volume  ratio  in  the  channel,  the  effective  diffusion 
coefficient  of  gas  and  liquid  phase  methanol  in  the  channel  can  be  described  as 


nMeOH 
UA-CH,l,eff  ' 


(1  -  a) 

*A 


nMeOH  nMeOH  _  a  nMeOH 
Ul  anQ  UA-CH,g,eff  ~  ug 


(36) 


Consequently,  the  two-phase  effective  diffusion  coefficient  in  the  mass-transfer  coefficient  as  shown  in  Eq.  (17),  is  properly 

described. 

As  a  result,  two  pressure  equilibrium  statements  at  the  anode  channel  and  the  diffusion  layer  are  set  up  with  undetermined  two-phase 
parameters,  S  and  a.  The  solution  can  be  found  from  knowledge  of  the  average  concentration  distribution  of  each  species,  namely 

and  C^%L  for  the  diffusion  layer,  and  and  C^°2H  for  the  anode  channel.  The  concentration  distribution  of  methanol  can  be  easily 

calculated  from  the  equations  derived  above  by  their  integration  over  the  entire  region.  The  derivation  of  the  concentration  distribution 
for  carbon  dioxide  will  be  presented  in  the  following  section. 


2.4.  Carbon  dioxide  mass  transport  at  anode 

Carbon  dioxide  is  a  product  of  the  fuel  cell  reaction.  Since  the  PEM  is  assumed  to  be  impermeable  to  the  gas,  there  is  no  internal  leaking 
flux  of  carbon  dioxide.  Therefore,  the  mass  flux  of  carbon  dioxide  travels  from  the  catalyst  layer  to  the  anode  channel  through  the  diffusion 
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layer  and  is  related  only  to  the  methanol  flux,  as  in  Eqs.  (14)  and  (23),  for  regions  with  or  without  under-rib  convection,  namely: 

(X)  (37) 

L 

Accordingly,  the  carbon  dioxide  concentration  in  the  anode  channel  £H[x)  can  be  found  from  the  mass  conservation: 


\jC02  _  atC02  _  1  j  _  1  r3  rMe0H  I 

'A—CH—DL  ~  " A-DL  ~  gpM  ~  Qp  1  A  '-A-DL  \y=tA_DL' 


tA-CHVA-CH 


dx 


-  —r3  c 

—  ^  ^  1  A  , 


■MeOH  I 


6 F  A  A~DL  b=tA-DL 


(x) 


(38) 


The  carbon  dioxide  concentration  at  the  interface  between  the  channel  and  the  diffusion  layer  CA°'foL  |  (x)  is  then  determined  by 

mass-transfer  transport: 


uC02  ( rC02  I  rD02  ^ ^  r> 3  /-Me 

nm,eff '  e A-DL  L=0  L A-CHJ  “  g p1  A  ^A-, 


-MeOH  I 
-DL  v_| 


y=tA-DL 


C X ) 


(39) 


while  the  concentration  in  the  channel/diffusion  layer  can  be  derived  as  in  previous  process  for  methanol  derivation  in  regions  with  or 
without  under-rib  mass  transport. 


2.4 A.  Under-rib  mass  transport  region 

The  constant  carbon  dioxide  flux  is  driven  by  a  concentration  gradient  as  well  as  a  pressure  difference  between  neighboring  channels, 
but  in  the  reverse  direction  to  the  methanol  flux,  i.e.: 


M C02  _  r»C02  dCA^)L  ,  r2cC02 

"A-DL  ~  UA-DL,ef(  dyA_DL  +  1  A  '-A-DL 


(40) 


where  the  effective  porosity  is  used  to  calculate  the  effective  diffusion  coefficient  of  carbon  dioxide,  D^°2)L  ^  =  [(1  -  S)£C]TcD^0 2.  The 
carbon  dioxide  distribution  in  the  diffusion  layer  can  be  derived  as 


CCA%L  = 


( ^ A-DL  | y= o  C A-DL  _  )exP( (  ^A  / DA-DL , eff  ^ A-DL  )  +  ^A-DL  \ y=tA_Di ,  ^ A-DL  | y=Q eXP( (  ^A  / DA-DL , eff  fa-DL ) 


y=tA-DL 


1  -  exp((rA2/  -  DcA%L  eff)tA-DL) 


(41) 


The  carbon  dioxide  flux  is  then  given  by 


NCa%l  = 


y=tA-DL 


-  n  CA%L\y=oe^-rA/DA%L,eff  W) 


\ -exp{(-r^/DcA%Leff)tA_DL) 


(42) 


2.4.2.  No  under-rib  mass  transport  region 

If  under-rib  convection  is  considered,  the  carbon  dioxide  flux  will  be  treated  as  a  diffusion  flow  linearly  distributed  across  the  diffusion 
layer,  i.e.: 


/«co2  _  Dco2  dCA?fa  _  Dco2 

*A—DL  -  u A-DL, eff  dyA_DL  ~  A-DL, eff 


CA-h\y=t,  D,  CA-bL\y=0 


tA-DL 


and, 


CCA%L  = 


^A-hly^u  CA-h 


y= o 


tA-DL 


yA-DL  +  CCA%L  |y=0 


(43) 


(44) 


Consequently,  the  average  concentrations  carbon  dioxide  in  the  anode  diffusion  layer  C^L  and  the  channel  C^H  are  now  available  for 
the  two-phase  model. 


2.5.  Oxygen  mass  transport  at  cathode 


The  situation  at  the  cathode  is  more  simple  without  the  two-phase  effect.  The  mass  balance  of  oxygen  at  the  cathode  side  is  similar  to 
that  of  methanol  at  the  anode  with  the  corresponding  parameters  in  their  gas  phase  value.  The  oxygen  supplied  at  the  cathode  catalyst 
site  combines  with  electrons  to  provide  electricity  generation,  and  also  partially  reacts  with  crossovered  methanol  and  thereby  reduces 
the  mixed  potential  of  the  cathode.  Thus 


N! 


1  3 

r°2  —  M°2  —  __  jV  _l  _  MJ 

C-CH-DL  ~  nC-DL  ~  4 pJC  +  9 


MeOH 


(45) 


where, 


DL  |  v=i 


y=k-DL 


rl  =  (j°}/C°})exp((acF/RT)ric) 


(46) 
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2.5 A.  Under-rib  mass  transport  region 

Oxygen  flux  is  driven  by  a  concentration  gradient  and  a  pressure  drop  between  adjacent  channels: 

NC-DL  =  ~DC-DL,eff  +  r^C°~DL 


32  Vc-ch^c-dl 


[2nlv  +  (2n-\)lh-2x] 


where  the  effective  diffusion  coefficient  is  D°fDL  =  (£C)TcD° fDV 

Following  the  same  derivation  of  anode  methanol  distribution,  the  oxygen  distribution  in  the  cathode  diffusion  layer  and  its  concentra¬ 
tion  at  the  catalyst  layer  can  be  expressed  as 


ce-DL  = 


LC-Dl|y=o  LC-DL 


)eXP( (  rl  C-DL , eff  ^ C~DL  )  +  ^ C-DL  y=tc_DL  ^C-DL  y=o eXP( (  ^ C  C-DL , eff  ) fC-DL  ) 


l-exp((r^/D°2DL^)tc_D1) 


C-DL  |  y=tc_DL 


_  ^w,U-rib  C02  \ 

~  ^C-DL  ^C-DL  |y=0  +  y 

where  the  equilibrium  constants  are 


rc  exp((rc1  /D°lDLeff)tC-DL) 


1  -  exp((rc’  /D°2_DLeff)tc-DL)  4F 


\  -  expd^  /D°2_DLeff)tc_DL)  4 F 

The  carbon  dioxide  concentration  at  the  interface  between  the  channel  and  the  diffusion  layer  can  be  derived  from  mass-transfer, 

N C-CH-DL  =  ^m,eff^C-CH  ~  ^C-DL  |y=o^ 

C°dDL  U  =  ^-CH^DL^CH  +  "m*™  (50) 

where  the  equilibrium  constants  are 

\i°2 

qjW.  U-rib  _  _ m,eff _ 

c-cn-DL  -  _  exp((rc,  /D02_dl  )tt_Dl)]/[i  _  exp((rc3 /D°2  )rc_Dl)]  +  Hh 

(51) 

^  = _ rc^7[i-exP((r(3/Dg2DLeff)tc_DL)] _ 

n  ^7-Dlrib  -  exp((rcVDg2DL  eff)tc_DL)]/[\  -  exp((ri/D°dDLeff)tc_DL)\  +  h%eff 

Finally,  the  oxygen  concentration  along  the  channel  C°fCH(x)  are  derived  by  considering  mass  conservation  and  distribution  in  the 
diffusion  layer,  i.e.: 

dc°  2 

-k-CHVC-CH-^  =  h%eff(  1  -  ^-CH-DL^C-CH  ~  ^ 

2.5.2.  No  under-rib  mass  transport  region 

Darcy-type  convection  flux  is  neglected  in  the  oxygen  flux.  Therefore,  the  oxygen  concentration  is  found  to  be  linearly  distributed  along 
the  diffusion  layer,  i.e.: 


N°2,  =  -D°2 


C-DL  ~  "C-DL, eff  dyc 


dC°-DL  =-  D°2, 


((-C-Di  ly=tc_D1  C°~DL  y=f 


Consequently, 


(CC-Dl|y=tc  DL  CC-DL  ) 

ro2  _  y  C~DL  y= o 

C~DL  tc-DL 


yc-DL  +  c°2dl  ^=q 


rO 2  I  _  xT/W/o, U-rib  ro2  I  _  NMeOH 

^C-DL\y=tc_DL  -  C-DL  '  ^C-DL\y=o  m 

r°2  I  -  I o,U—rib  ro2  ,  xjp'  jsjMeOH 
u C-DL  |y=o  ~~  C-CH-DL  '  u C-CH  ^  *  '  iym 


Table  1 

Physicochemical  properties. 


Y.  Yang,  Y.C.  Liang  /  Journal  of  Power  Sources  194  (2009)  712-729 


721 


Parameters 

Symbol 

Value 

Anode  inlet  pressure  (Pa) 

P Ain 

1.01325  x  105 

Cathode  inlet  pressure  (Pa) 

Pan 

1.01325  x  105 

Oxygen  input  concentration  at  cathode  (mol cm-3) 

Cc 

in 

0.21  x  10-6  { PcmIRT ) 

Permeability  of  PEM  (cm2) 

K 

1.17  x  1019  exp(-19098/T) 

Permeability  of  anode  diffusion  layer  (cm2)  [31] 

k A-DL 

1  x  19-9 

Permeability  of  anode  diffusion  layer  (cm2)  [31] 

K-ol 

1  x  19-9 

Electro-osmosis  drag  of  water  (cm3  mol-1 )  [35] 

nE$S 

2.5 

Electro-osmosis  drag  of  methanol  (cm3  mol-1)  [35] 

yMeOH 

^eod 

45 

Viscosity  of  liquid  water  (Pa  s)  [39] 

Va 

0.458509  -  5.30474  x  10-3T+ 2.31231  x  10-5T2  -  4.49161  x  10-8T3 +3.27681  x  10 ~nP 

Viscosity  of  gas  (Pa  s)  [39] 

[XC 

2.03  x  10-5 

Diffusion  coefficient  of  liquid  methanol  (cm2  s-1)  [40] 

£jMeOH 

(10(-5.4163-(999. 778/7)))  x  104 

Diffusion  coefficient  of  gas  methanol  (cm2  s-1 )  [40] 

£)MeOH 

-6.954  x  10-2  +  4.5986  x  10-47+  9.4979  x  10-7 7 2 

Diffusion  coefficient  of  carbon  dioxide  (cm2  s-1 )  [18] 

frq  n 
O 

0.3 

Diffusion  coefficient  of  oxygen  (cm2  s-1 )  [41  ] 

D°2 

1.77(7/273)1823(  10.13/Pc, in) 

Porosity  of  PEM 

0.3 

Porosity  of  anode  diffusion  layer 

£a 

0.6 

Porosity  of  cathode  diffusion  layer 

sc 

0.7 

Tortuosity  of  PEM 

Tm 

1.8 

Tortuosity  of  anode  diffusion  layer 

r a 

1.5 

Tortuosity  of  cathode  diffusion  layer 

r  c 

1.5 

Proton  conductivity  of  PEM  (S  cm-1 )  [42] 

Km 

0.123 

Reference  anode  exchange  current  density  (A cm-3)  [18] 

,-MeOH 

'ref 

9.425  x  10-3  exp((35570)/R((  1  /353 )-( 1/7))) 

Reference  cathode  exchange  current  density  (A cm-3)  [18] 

j°  2 
'ref 

4.222  x  10-6  exp((73200/R)((l/353)  -  (1/7))) 

Reference  concentration  of  methanol  (mol  cm-3) 

r  MeOH 
ref 

0.001 

Reference  concentration  of  oxygen  (mol  cm-3) 

C°  2 

ref 

0.002 

Anodic  transfer  coefficient 

OLa 

0.54 

Cathodic  transfer  coefficient 

Oic 

0.76 

Sherwood  number  [18] 

Sh 

2.693 

Interfacial  tension  (S cm-1 )  [37] 

o 

6.25 

Contact  angle  (°) 

0 

110 

Contact  resistance  (C2  cm2) 

R contact 

0.01 

where  the  equilibrium  constants  are 


tfjW / o  ,U -rib 
C-DL 


DC— DL,  eff  /^C—DL 


D8W/tc-“+(rc/4F) 


and  ^ 


'w/o,U-rib 


3/2 


C-DL 


w/o,U-rib 


h°  2 
m,eff 


C-CH-DL 


and 


w/o,U-rib 


DclDLeff/tc-DL  +  ^/A-F) 

(1/4 F)r^7l°'LU~rib  -  (3/2) 


-tc-CHVC-CH 


0  /4F)r^!°DLu-rib  +  h°m2eff  C~CH~DL  (i  /4F)r^!°DLu-rib  +  h°m2e 

the  cathc 

C-CH  _h02  M  _lr/w/o,U-rib)co2  _h02  tf/'KiMeOH 
~  nm,eff'1  * C-CH-DL  C-CH  nm,eff  * 


(55) 


(56) 


lC^C-DL  K L  /  -*1  c  ^  C—DL  ^  ,lm,eff 

Similarly,  the  oxygen  concentration  in  the  cathode  channel  and  its  distribution  in  the  diffusion  layer  can  then  be  determined,  namely: 

dC°  2r 


dx  m,eff '  C-CH-DL  '  C-CH  m,eff 

The  cell  current  density  is  then  calculated  by  integration  of  the  local  current  density  at  cathode. 

f^total 

Icell  — 


rLtotal 

■■  /  rc2c°/Di(x)dx 

JO 


(57) 


(58) 


2.6.  Cell  voltage 


Once  the  current  density  and  the  overpotentials  of  the  anode  and  the  cathode  are  obtained,  the  cell  voltage  can  be  determined  by 
subtracting  the  activation  losses  at  the  two  electrodes,  as  well  as  the  ohmic  contact  potential  and  PEM  potential  from  the  thermodynamic 
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Fig.  2.  Validation  of  present  model  against  experimental  data  for  a  5  cm2  single-cell  DMFC:  (a)  feeding  with  different  concentration  methanol  solutions,  0.125,  0.25,  0.5,  and 
1  M;  (b)  under  different  operating  temperatures,  25, 40  and  60  °C,  with  1  M  feeding  methanol  solution.  Flow  rate  of  methanol  in  anode  channel  is  2  ml  min-1.  Air  is  supplied 
at  cathode  with  flow  rate  of  300  ml  min-1  in  all  cases. 


equilibrium  potential.  Hence 

Vcell  =  V°2  -  VrH  -VA-VC-  hell  ^  ~  kellRcontact  (59) 

Km 

where  and  v^e0H  are  the  thermodynamic  equilibrium  potentials  of  oxygen  reduction  and  methanol  oxidation,  their  difference  is  set 
as  1.21  V  for  the  theoretical  thermodynamic  cell  voltage  adopted  in  this  study. 

3.  Experimental  setup 

3  A.  DMFC 

A  single-cell  DMFC  that  consisted  of  a  membrane  electrode  assembly  (MEA)  sandwiched  between  two  graphite  blocks  (POCO  Graphite 
Inc.,  Decatur,  TX)  was  used  for  performance  measurement.  The  active  area  of  the  MEA  was  5  cm2  and  a  Nation®  117  membrane  was 
employed.  The  anode  catalyst  loading  on  the  membrane  was  4  mg  cm-2  of  Pt-Ru  black  with  an  atomic  ratio  of  1:1.  At  the  cathode  side, 
only  Pt  was  applied  with  the  same  loading  of  4  mg  cm-2.  Before  assembly  into  a  single-cell,  the  MEA  was  pretreated  to  improve  its  water 
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Fig.  3.  Modelling  results  of  methanol  concentration  distribution  within  anode  channel,  diffusion  layer  and  PEM  for  different  input  methanol  concentrations,  0.125,  0.25,  0.5, 
and  1  M,  at  Vceu  =  0.35  V  and  40  °C,  2  ml  min-1  for  MeOH,  300  ml  min-1  for  air. 


content  by  rinsing  in  50  °C  deionized  (DI)  water  for  5  h  and  then  storing  overnight  in  DI  water  at  room  temperature.  The  flow-fields  for 
fuel  transport  and  distribution  were  fabricated  from  graphite  blocks  with  a  SSC  pattern.  The  channel  had  a  square-shaped  cross-section. 
The  entire  flow-field  consisted  of  14  parallel  channels.  The  dimensions  of  each  part  in  the  flow-field  are  described  in  Fig.  1.  With  two 
polytetrafluoroethylene  (PTFE)  coated  fiberglass  gaskets,  the  flow-field  plates  were  fastened  to  the  MEA,  with  carbon  cloth  (E-TEK  B-l,  Type 
A)  and  E-TEK  ELAT®  (LT1400-W)  as  the  diffusion  layers  at  the  anode  and  cathode  respectively.  The  current-collectors  on  both  sides  were 
gold  plated  to  reduce  electrical  resistance. 

3.2.  Measurement  conditions 

After  assembly,  the  DMFC  hardware  was  placed  in  a  temperature-controlled  chamber  (ETAC,  HISPEC  FIT210)  for  performance  measure¬ 
ments.  The  temperature  was  monitored  with  a  thermometer  (FLUKE  179)  attached  to  the  flow-field  plates.  At  the  anode  side,  the  aqueous 
methanol  solution  was  supplied  by  a  digital  peristaltic  pump  (Longer,  LEAD-1)  with  an  adjustable  flow  rate.  At  the  same  time,  the  oxy¬ 
gen  was  supplied  by  a  diaphragm  air-pump  (FIAILEA,  AC0-9610)  from  the  ambient  environment.  An  electronic  load  (Chroma  63106)  was 
employed  to  record  the  current-voltage  (/-V)  curves.  Before  each  measurement,  the  DMFC  was  first  operated  under  a  high  cell  voltage 
output  (0.5  V)  and  a  high  current  output  (0.9  A)  at  60  °C  for  5  h  to  activate  fully  the  MEA. 

4.  Results  and  discussion 

4.1.  Model  validation 

The  experimental  data  for  model  validation  was  collected  with  the  above  setup  under  different  operating  conditions,  i.e.,  variable 
methanol  concentrations  (0.125,  0.25,  0.5,  and  1  M)  and  operating  temperatures  (25,  40,  and  60).  The  flow  rate  of  methanol  is  set  at 
2  ml  min-1 ,  and  air  flow  rate  is  300  ml  min-1 .  The  physical  parameters  used  in  the  proposed  model  are  tabulated  in  Table  1 .  In  Fig.  2a,  the 
modelling  results  of  the  cell  l-V  curves  at  different  methanol  concentrations  and  40  °C  operating  temperature  are  verified  and  compared 
with  experimental  data.  The  results  show  a  reasonable  agreement  between  the  two  sets  of  data.  The  same  is  also  observed  in  Fig.  2b  for 
DMFC  performance  at  different  operating  temperatures  with  1  M  methanol  solution.  Thus,  the  mathematical  model  developed  in  this  study 
is  believed  to  be  realistic  for  DMFC  performance  analysis  and  prediction. 

For  various  methanol  concentrations  at  40  °C  and  a  cell  output  voltage  of  0.35  V,  the  corresponding  average  methanol  concentration 
within  the  channel,  diffusion  layer  and  PEM  are  profiled  in  Fig.  3.  The  distribution  of  the  methanol  concentration  across  the  anode  diffusion 
layer  and  PEM  are  also  depicted  in  the  expanded  figure.  The  non-linear  methanol  distribution  in  the  diffusion  layer  implies  that  under¬ 
rib  mass  transport  improves  mass  transport  in  this  layer.  The  relatively  low  methanol  concentration  in  the  catalyst  layer  indicates  that 
the  electrochemical  reaction  at  the  anode  is  still  restricted  by  methanol  transport,  even  though  excess  methanol  is  supplied  in  the  anode 
channel.  In  all  regions,  the  diffusion  layer  plays  a  major  role  in  the  total  mass-transfer  resistance.  Higher  methanol  concentration  at  the  input 
enhances  mass  transport.  On  the  other  hand,  it  also  induces  higher  methanol  crossover  flux  in  the  PEM,  as  seen  in  Fig.  3.  The  limiting  mass- 
transport  phenomenon  occurs  earlier  with  lower  methanol  concentration  at  the  input,  as  indicated  in  Fig.  4.  The  methanol  concentration 
in  the  catalyst  layer  is  nearly  zero  with  0.125  M  of  input  concentration.  Consequently,  the  crossover  flux  in  the  PEM  is  much  lower  since 
most  of  the  methanol  is  consumed  for  the  electricity  generating  reaction  on  the  catalyst. 
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Fig.  4.  Leaking  current  density  caused  by  methanol  crossover  for  different  input  methanol  concentrations,  0.125,  0.25,  0.5,  and  1  M,  with  2  ml  min-1  flow  rate  at  40  °C. 


4.2.  Analysis  of  methanol  crossover 

Methanol  crossover  is  driven  by  diffusion,  electro-osmosis  and  convection.  These  three  factors  contribute  differently  under  different 
operating  conditions.  Nevertheless,  all  three  depend  on  the  local  methanol  concentration  in  the  PEM.  Hence,  the  maximum  crossover  flux  is 
expected  at  the  open-circuit  voltage  where  the  highest  mass-transfer  of  un-reacted  methanol  occurs.  When  the  current  density  increases, 
the  excess  methanol  at  the  interface  between  the  anode  catalyst  layer  and  the  PEM  decreases.  Accordingly,  methanol  crossover  reduces  to 
almost  zero  after  the  current  density  approaches  its  limiting  point  where  nearly  all  methanol  is  consumed  by  the  electrochemical  reaction. 
This  is  depicted  in  Fig.  4  by  converting  the  methanol  crossover  flux  into  the  induced  leakage  current  density,  i.e.: 

Ileak  =  6  FN™e0H  (60) 

The  high  leakage  current  density,  which  occurs  in  the  low  output  current  density  region,  brings  more  severe  voltage  reduction,  as 
shown  in  Fig.  2a.  The  two  partial  fluxes,  namely,  diffuse  and  electro-osmosis  drag,  in  methanol  crossover  have  been  investigated  and 
are  shown  in  Fig.  5  for  0.5  and  1  M  input  methanol  densities,  respectively.  Due  to  the  lack  of  a  significant  pressure  difference  between 
the  anode  and  the  cathode,  the  convection  contribution  is  neglected.  The  result  indicates  that  diffusion  mainly  dominates  in  methanol 
crossover  except  in  the  region  near  the  limiting  current  density.  In  this  region,  the  methanol  concentration  at  the  interface  of  the  catalyst 
layer  and  the  PEM  decreases  due  to  the  mass  transport  limitation  and  high  current  output.  Therefore,  the  diffuse  methanol  flux  governed 
by  the  concentration  gradient  in  the  PEM  is  markedly  suppressed.  The  electro-osmosis  drag  flux,  on  the  other  hand,  does  not  decrease 
significantly  because  it  is  proportional  not  only  to  the  decrease  in  methanol  concentration  but  also  to  the  increase  in  cell  current  density. 
Once  the  current  density  reaches  its  limiting  point,  the  diffuse  flux  becomes  zero  since  there  is  no  concentration  gradient.  The  electro¬ 
osmosis  drag  flux  of  methanol  also  reduces  to  zero  because  there  is  no  methanol  content  in  the  electro-osmosis  water  flux  transferring 
through  the  PEM. 

The  temperature  effect  on  the  methanol  concentration  induced  crossover  is  presented  in  Fig.  6.  Methanol  crossover  gives  us  to  a  voltage 
reduction  at  low  cell  current  density.  The  phenomenon  is  enhanced  by  raising  the  operating  temperature  from  40  to  60  °C.  It  is  thought  that 
temperature-enhanced  methanol  transport  brings  a  relatively  higher  methanol  concentration  at  the  interface  between  the  catalyst  layer 


Fig.  5.  Comparison  between  diffusion  and  electro-osmosis  drag  partial  flux  in  methanol  crossover  flux  for  0.5  and  1  M  input  methanol  concentrations,  with  2  ml  min-1  flow 
rate  at  40  °C. 
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Fig.  6.  Cell  voltage  loss  caused  by  methanol  crossover  in  low  current  density  regions  with  methanol  concentration  increase  (0.5-1  M)  with  2  ml  min-1  flow  rate  under  different 
operating  temperatures,  (40  and  60  °C). 


and  PEM.  At  higher  operating  temperatures,  a  higher  diffusion  coefficient  applies  to  both  liquid  and  gas  phase  methanol.  At  the  same  time, 
the  saturated  methanol  concentration  in  the  gas  phase  also  increases  to  give  fast  methanol  transport. 

4.3.  Effects  of  two-phase  and  under-rib  mass  transport 

To  investigate  the  two-phase  mass  transport  effect  on  cell  performance,  the  I-V  curve  for  1  M  input  methanol  at  40  °C  was  simulated. 
The  liquid  saturation  S  is  assumed  to  be  unity  and  thereby  represents  only  liquid  phase  transport  at  the  anode.  The  result  is  shown  in 
Fig.  7  for  the  case  of  an  under-rib  mass  transport  effect.  The  limiting  current  density  is  147.0  mA  cm-2  which  is  much  smaller  than  that 
for  two-phase  mass  transport,  viz.,  200.6  mAcm-2.  Thus  reveals  that  methanol  transport  in  the  gas  phase  is  important  for  a  sufficient  fuel 
supply  to  the  reaction  sites.  Moreover,  the  cell  voltage  reduction  induced  by  methanol  crossover  is  noted  at  a  low  current  density  where 
the  liquid  phase  transport  situation  gives  a  higher  cell  output  voltage  compared  with  that  for  two-phase  transport.  The  phenomenon  is 
attributed  to  enhanced  methanol  transport  in  the  gas  phase  that  in  more  severe  methanol  crossover.  Accordingly,  the  theoretical  limiting 
current  density  can  be  calculated  by  setting  the  methanol  concentration  to  zero  at  the  interface  between  the  anode  catalyst  layer  and  the 


Fig.  7.  Influence  of  two-phase  flow  phenomenon  on  DMFC  performance,  with  and  without  under-rib  mass  transport  consideration  effects,  1  M  input  methanol  solution  with 
2  ml  min-1  flow  rate  at  40  °C  operating  temperature. 
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PEM,  i.e.: 


^lim  — 


6  F 

(i/Ce0H)  +  (T77^j[T  -  (i /exP((ry2/DHf  )t^_DL))] 


rMeOH 

^A-CH 


(61) 


With  an  input  methanol  concentration  of  1  M,  the  calculated  limiting  current  density  is  153.7  mAcm-2,  which  is  consistent  with  the 
modelling  result  given  in  Fig.  7. 

By  further  ignoring  the  under-rib  mass  transport  effect,  the  above  equation  can  be  simplified  as 


6  F 

(l/h^0H)  +  (^_DL/Df^H) 


rMeOH 
V A-CH 


(62) 


Consequently,  the  limiting  current  density  has  a  relatively  smaller  value  of  135.5  mA  cm-2  with  both  effects  ignored,  which  is  also 
reasonably  close  to  the  modelling  result  of  138.4  mA cm-2.  The  difference  with  and  without  under-rib  mass  transport  effect,  implies  that 
the  facilitated  methanol  transport  toward  the  catalyst  layer  is  influenced  by  under-rib  mass  transport.  The  significance  of  under-rib  mass 
transport  on  performance  is  however,  less  than  that  caused  by  the  two-phase  effect.  Performance  data  with  the  under-rib  mass  transport 
effect  subtracted  from  the  two-phase  mass  transport  are  presented  in  Fig.  7. 
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Fig.  8.  Flow-field  structure  and  channel  cross-section  effects  on  DMFC  limiting  current  density  and  maximum  power  output:  (a)  under-rib  width  increases  with  same  cross- 
section  area;  (b)  under-rib  width  increases  with  same  channel  depth.  Entire  flow-field  area  and  total  parallel  channel  number  unchanged  in  both  cases.  1  M  input  methanol 
solution  with  2  ml  min-1  flow  rate  at  40  °C  operating  temperature. 
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Fig.  9.  Influence  of  operating  parameters  on  DMFC  limiting  current  density  and  maximum  power  output  (a)  input  methanol  concentration  increases  from  0.25  up  to  6  M;  (b) 
operating  temperature  increases  from  25  to  60  °C.  Flow  rate  of  methanol  solution  varies  from  1  to  30  ml  min-1  for  both  cases. 


4.4.  Flow-field  structure 

As  mentioned  earlier,  the  diffusion  layer  imposes  a  significant  effect  on  the  mass  transport  resistance  in  a  DMFC.  An  investigation  of  its 
influence  on  cell  limiting  current  density  and  maximum  power  density  was  carried  out  by  changing  the  under-rib  width  and  cross-section 
while  keeping  the  same  total  active  area  and  channel  parallel  number  in  the  field.  Fig.  8a  depicts  the  results  obtained  from  the  modified 
flow-field  structure.  In  the  modification,  the  under-rib  width  between  two  adjacent  channels  was  increased  from  0.01  to  0.14  cm,  while 
the  cross-sectional  area  was  kept  constant  at  0.075  cm  x  0.075  cm.  With  the  total  flow-field  area  and  channel  parallel  number  unchanged, 
the  contact  length  between  the  channel  and  diffusion  layer  then  decreases  from  0.1105  to  0.0455  cm.  At  the  same  time,  the  cross-sectional 
shape  of  each  channel  will  also  change.  The  change  in  cross-section  of  two  adjacent  channels  with  a  connecting  under- rib  is  given  in  Fig.  8a 
for  a  horizontal  shape,  square  shape,  and  a  vertical  shape.  Both  the  limiting  current  density  and  maximum  power  density  reduce  with 
the  increasing  under-rib  width  due  to  increased  mass  transport  resistance.  Since  the  mass-transfer  at  the  channel|diffusion  layer  interface 
depends  on  the  contact  area  and  the  diffusion  flow  mainly  occurs  underneath  this  open  space,  it  is  concluded  that  the  performance 
degradation  is  primarily  caused  by  inefficient  mass-transfer  and  reduced  diffuse  flux.  Even  though  the  flow  velocity  in  the  channel  does 
not  change  for  constant  cross-sectional  area,  the  pressure  difference  between  adjacent  channels  becomes  higher  with  decreasing  channel 
hydraulic  number.  The  under-rib  mass  transport  will  then  be  enhanced  by  the  increasing  pressure.  The  contrasting  contributions  on  cell 
performance  from  mass-transfer/diffusion  and  under-rib  convection  reveal  that  effective  mass-transfer  from  the  channel  and  its  diffusion 
in  the  diffusion  layer  are  dominating  factors  the  optimum  mass  transport  for  electrochemical  reaction.  Given  in  the  steep  drop  in  limiting 
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current  density,  the  drop  in  maximum  power  density  in  Fig.  8a  is  the  result  of  heavy  methanol  crossover  in  the  low  mass-transport  resistance 
region. 

Another  analysis  is  given  in  Fig.  8b  with  a  similar  structure  modification  as  that  in  Fig.  8a  while  maintaining  the  channel  height  at 
0.075  cm,  as  adopted  in  the  experimental  hardware,  where  the  flow  velocity  in  the  channel  increases  with  increasing  under- rib  width.  The 
higher  flow  velocity  results  in  more  effective  mass-transfer  at  the  channel] diffusion  layer  interface.  It  also  facilitates  mass  transport  in  the 
under-rib  region  by  improving  the  pressure  difference  between  adjacent  channels.  Therefore,  the  performance  of  the  DMFC  is  promoted 
with  a  wider  under-rib  width.  Nevertheless,  the  improvement  is  moderate  as  the  contact  area  between  the  channel  and  diffusion  layer 
shrinks  quickly  and  thus  limits  the  effective  mass-transfer. 

Both  analyses  show  the  influence  of  the  open  ratio  in  the  flow-field  channel  regions  for  effective  mass-transfer  at  the  channel] diffusion 
layer  interface  and  its  diffusion  flow.  The  reduced  contact  between  the  flow-field  rib  and  diffusion  layer  also  brings  a  high  electrical 
resistance  in  electron  transport.  As  a  result,  a  trade-off  should  be  considered  to  achieve  an  optimum  flow-field  design. 

4.5.  Influence  of  operating  parameters 

Variations  in  the  limiting  current  density  and  maximum  power  density  as  functions  of  methanol  concentration  and  temperature  for 
different  methanol  flow  rates  are  presented  in  Fig.  9a  and  b,  respectively.  For  all  cases,  the  limiting  current  densities  increase  with  methanol 
flow  rate  at  the  anode.  There  is,  however,  no  further,  improvement  at  flow  rates  above  5  ml  min-1.  A  similar  trend  can  also  be  observed 
for  the  power  density,  which  changes  moderately  compared  with  that  in  the  limiting  current  density.  Under  a  slow  flow  rate,  the  low 
liquid  velocity  in  the  anode  channel  leads  to  a  lower  mass-transfer  at  the  interface  between  the  channel  and  diffusion  layer.  The  pressure 
drop  across  the  SSC  is  also  smaller,  and  thus  induces  a  weaker  pressure  that  drives  under-rib  convection.  Consequently,  the  under-rib 
mass  transport  is  suppressed.  As  a  result,  the  diffusion  layer  exhibits  a  higher-mass-transfer  resistance  against  the  methanol  supply  at  the 
reaction  sites.  Although  flow  rate  adjustment  will  improve  the  situation,  a  large  increase  in  flow  rate  will  bring  no  significant  improvement 
in  methanol  transport  due  to  the  absence  of  sufficient  fuel  residence  time. 

The  limiting  current  density  is  found  to  increase  with  higher  methanol  concentration.  On  the  other  hand,  the  increase  in  maximum 
power  density  becomes  less  remarkable  after  2  M  and  drops  at  6  M  in  Fig.  9a.  The  results  indicate  that  the  enhanced  mass-transfer  for 
electrochemical  reactions  due  to  increased  input  methanol  concentration  also  makes  the  methanol  crossover  more  severe  and  eventually 
results  in  a  poor  power  density.  By  contrast,  an  increase  in  operating  temperature  facilitates  an  improvement  in  limiting  current  density  and 
power  density.  Effective  methanol  transport  is  ensured  under  a  high  operating  temperature  because  of  its  high  diffusion  coefficient.  At  the 
same  time,  the  reaction  kinetics  are  also  improved,  and  the  methanol  crossover,  which  results  from  excess  methanol  without  reaction,  can 
be  reduced.  Overall,  the  control,  of  input  methanol  concentration  has  a  more  significant  effect  than  the  suppression  of  methanol  crossover 
in  achieving  an  improvement  in  DMFC  performance. 

5.  Conclusion 

In  this  study,  a  complexity-reduced  mathematical  model  has  been  developed  for  a  liquid-feed  DMFC.  The  model  involves  two  significant 
mass-transport  phenomena,  namely:  (i)  the  two-phase  transport  effect  due  to  carbon  dioxide  generation  at  the  anode,  and  (ii)  the  under¬ 
rib  mass  convection  effect  in  the  presence  of  a  pressure  drop  along  a  single  serpentine  fuel  channel.  This  model  provides  a  more  realistic 
performance  analysis  of  DMFC  design.  The  results  obtained  from  the  model  are  in  good  agreement  with  experimental  data.  Thus,  proposed 
model  serves  as  a  useful  tool  for  the  optimization  of  DMFC  design  and  performance. 
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